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ABSTRACT 



There has been increasing evidence that at least some GRBs are emission 
beamed. The beamed GRB afterglow evolution has been discussed by several 
authors in the ultra-relativistic case. It has been shown that the dynamics 
of the blast wave will be significantly modified by the sideways expansion, 
and there may be a sharp break in the afterglow light curves under certain 
circumstances. However, this is true only when the fireball is still relativistic. 
Here we present an analytical approach to the evolution of the beamed GRB 
blast wave expanding in the surrounding medium (density n oc r~ s ) in the 
non-relativistic case, our purpose is to explore whether the sideways expansion 
will strongly affect the blast wave evolution as in the relativistic case. We 
find that the blast wave evolution is strongly dependent on the speed of the 
sideways expansion. If it expands with the sound speed, then the jet angle 9 
increases with time as 9 oc Int, which means that the sideways expansion has 
little effect on the afterglow light curves, the flux F oc t s for s = and 

7a + l 

F oc t ~ for s = 2. It is clear that the light curve of s — 2 is not always 
steeper than that of s = 0, as in the relativistic case. We also show that if the 
expansion speed is a constant, then the jet angle 9 oc t, and the radius r oc t , 
in this case the sideways expansion has the most significant effect on the blast 
wave evolution, the flux F oc £ _ ( 5a_1 ) independent of s, and we expect that 
there should be a smooth and gradual break in the light curve. 

Key words: gamma-rays: bursts 
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1 INTRODUCTION 

The fireball model of 7-ray bursts led to the prediction of the afterglow emission that might 
be expected when the energetic shock wave encountered the surrounding medium. The sub- 
sequent X-ray and optical observations of the afterglow appeared to confirm the prediction 
of the simplest afterglow model (Wijers, Rees k Meszaros 1997). This model involved syn- 
chrotron emission from electrons accelerated to a power law energy spectrum in a relativistic 
blast wave. However, both geometry and environment can affect the evolution of 7-ray bursts 
afterglows (Panaitescu, Meszaros k Rees 1998). There has been increasing evidence that at 
least some GRBs are emission beamed or the surrounding medium is not uniform. A class of 
GRBs whose afterglow exhibited steeper than the normal power law decays (f v oc t~ 2 ) can 
well be explained by the jet -like geometry of the relativistic shock (Rhoads 1997,1999; Sari, 
Piran k Halpern 1999; Wei k Lu 1999), or the inhomogeneous surrounding medium models 
(especially the wind model, Dai k Lu 1998; Chevalier k Li 1999; Li k Chevalier 1999). The 
jet model is also supported by the steepening of the optical and radio light curves seen in 
GRB 990510 (Stanek et al. 1999; Harrison et al. 1999). 

The dynamical evolution of GRB fireballs and the emission features have been studied 
by many authors (e.g. Sari 1997; Meszaros, Rees k Wijers 1998; Wei k Lu 1998a, 1998b; 
Sari, Piran k Narayan 1998), but most of them considered the fireball being isotropic. The 
evolution of the beamed blast wave has been first discussed by Rhoads (1997, 1999). He 
has shown that the lateral expansion of the relativistic plasma causes that at some moment 
the surface of the blast wave starts to increase faster than due to the cone-outflow alone, 
then the blast wave begins to decelerate faster than without the sideways expansion since 
more intersteller medium has been swept up by the blast wave. He claimed that this effect 
will produce a sharp break in the GRB afterglow light curves. More detailed calculation, 
both numerical and analytical studies (Moderski, Sikora k Bulik 1999; Wei k Lu 1999) have 
found that unless the opening angle is very small or that lateral expansion is unimportant, 
a smooth and gradual transition is expected. However, these results are valid only when the 
blast wave is relativistic. As shown by Rhoads (1999) and Panaitescu k Meszaros (1999), 
the Lorentz factor T b at the radius r b where the sideways expansion becomes important is 
r 6 ~ ^775^0 where 9 is the initial half opening angle of the jet. To keep T b ^> 1 requires 
#0 < 0.1, so if the jet angle is not too small, the blast wave has become non-relativistic when 
sideways expansion is important, and the previous results are invalid. 
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Frail et al. (1999) reported on the results of an extensive observation of the radio after- 
glow of GRB970508, lasting 450 days after the burst. They have shown that the spectral 
and temporal radio behavior indicate that the fireball has undergone a transition to non- 
relativistic expansion at t ~ 100 days, and they find that the fireball may be initially a wide 
angle jet of opening angle ~ 30°. Therefore it is very interesting and important to study the 
behavior of the beamed blast wave in the non-relativistic case. Recently Huang et al.(1999) 
have calculated the evolution of jetted GRB ejecta numerically. 

Here we present an analytical approach to the evolution of the beamed GRB afterglow 
in the non-relativistic case, including both homogeneous medium and wind-shaped medium 
(the density n oc r -2 ), the main purpose is to explore whether the sideways expansion will 
strongly modify the blast wave evolution and the afterglow light curve behavior as in the 
relativistic case. In next section we consider the dynamical evolution of blast wave, in section 
3 we calculate the jet emission and afterglow light curve analytically, and finally we give 
some discussions and conclusions. 



2 DYNAMICAL EVOLUTION OF THE JET 

Now we consider the evolution of an adiabatic blast wave expanding in surrounding medium. 
Assuming that the medium density n oc r~ s , s = corresponds to the homogeneous medium, 
and s = 2 corresponds to the wind-shaped medium. For completeness, here we first outline 
the results of evolution in the relativistic case. 

2.1 relativistic case 

For energy conservation, the evolution equation of blast wave is 

T 2 N = const (1) 

where T is the bulk Lorentz factor, and N is the total baryon numbers swept up by the 
blast wave, N oc r 3 ~ s (l — cosOj) oc r 3_s 6> 2 for 9j <C 1, and 9j = 9 + 9', where 9 is the 
initial jet opening half-angle, 9' describes the lateral expansion, which can be simply written 
as 9' ~ Vgtco/ct, v s is the expanding velocity of ejecta material in its comoving frame, and 
t (t co ) is the time measured in the burster frame (comoving frame). For relativistic expanding 
material it is appropriate to take v s to be the sound speed v s = c/ Since the jet expands 
relativistically, there is the relation T oc r/T 2 , where T is the time measured in the observer 
frame. According to Wei & Lu (1999), for s = 0, we obtain 
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* \ T- 3 / 8 [l + (J)V2]-i/4 ; if t>T 5 
where T fe is the time at that moment the sideways expansion is important. Similarly, for 
s = 2, we have 

f r- 1 /4[l + (T ) l/4 ] -l/2 > ifT<Tfc 
^ I T- 1 /4[l + (J)l/2]-l/2 ? if T>Tfo 

so we see that, for T < T 6 , T oc T~ 3 / 8 and r oc T 1 / 4 for s = 0, and T oc T" 1 / 4 and r oc T 1 / 2 
for s = 2, while for T > T 5 , both s = and s = 2 give T oc T" 1 / 2 and r oc T°. 



2.2 non-relativistic case 

The evolution of the jetted blast wave in the non-relativistic case is 

f3 2 N = const (4) 

where (3 is the velocity of the blast wave in units of light speed c. Thus, for ^ < 1, we have 
[3 2 r 3 ~ s 6j = const. In this case the jet angle 6j can be written as 9j = 6o + I j^|> where (3 S 
is the spread velocity of ejecta material in units of light speed c, which cannot be simply 
determined. Here we consider several situations. 

(i) P s =sound speed Kirk & Duffy (1999) have derived the sound speed in the fluid 
2 (7-l)p 

where P is the pressure, p is the mass density, and 7 is the adiabatic index. According to 
Huang et al. (1999), the sound speed can be written as 

(6) 



2 7(7-l)(7-l) 



1+7(7-1) 

In the non-relativistic limit (7 ~ 1, 7 ~ 5/3), one gets c s = a/5/3/3, then we obtain the jet 
angle 6j = + ^ n (^)' wnere ^0 ~ ijvfl, t^R is the time when the blast wave turns from 
relativistic to non-relativistic. Therefore we see that since the jet angle 0j increases with 
time as the log relation, it will not strongly affect the evolution of the jet, the variation of 
the jet velocity is about (3 oc t~( 3 ^ s )/( 5 ~ s ) ; i.e., (3 oc t~ 3 / 5 for s = 0, and (3 oc t~ 4 / 3 for s=2. 

(ii) (3 S = const If the jet spreads at a constant speed, then for t ^> t NR , the angle 9j ~ 
^0 + 3^tj, and the evolution of (3 satisfies (3 5 ' s (9 + f) 2 oc r (3 ~ s) , so when (3 < ^(3 S 9 \ 
we get (3 oc r ~ const independent of s. We see that in this case it is similar to the 
situation that sideways expansion is important in the relativistic case. 
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(iii) In general, we assume that the ratio of the jet spread velocity and the blast wave 
velocity varies with time, (3 s /(3 = at q , then we obtain when t ^> t^R, the jet angle 9j ~ 
9 + q (^l s ) t q , and for t » ((5 - s) q0 /2a) 1/<? , the evolution of (3 satisfies (3 oc t~ ' , so 
(3 oc t"^^ 2 for s = and (3 oc t ^r 1 for s = 2. We see that q — 1 corresponds to the case (ii). 
Therefore it seems that the evolution of the blast wave strongly depend on the parameter q. 



3 THE EMISSION FROM THE JET 

Now we calculate the emission flux from the non-relativistic jet. Here we adopt the formu- 
lation and notations of Mao & Yi (1994). In our model the ejecta is flowing outwards in a 
cone with opening half angle 9j. For simplicity, we assume that the radiation is isotropic in 
the comoving frame of the ejecta and has no dependence on the angular positions within the 
cone. The radiation cone is uniquely defined by the angular spherical coordinates (9, 0), here 
9 is the angle between the line of sight (along z-axis) and the symmetry axis, and is the 
azimuthal angle. Because of cylindrical symmetry, we can assume that the symmetry axis of 
the cone is in the y — z plane. In order to see more clearly, let us establish an auxiliary coordi- 
nate system (x', y', z') with the z'-axis along the symmetry axis of the cone and the x' parallel 
the x-axis. Then the position within the cone is specified by its angular spherical coordinates 
9' and <fi' (0 < 9' < 9j, < 0' < 2ir). It can be shown that the angle 9 between a direction 
(9', (f)') within the cone, and the line of sight satisfies cosQ = cos9cos9' — sin9sin9' siri(p' . 
Then the observed flux is 

F(i/,9) = r d(j)' f 3 sin9'd9'cosQD z I'(vD- l ) r ^ (7) 
jo jo d 2 

where D = [T(l - pcosQ)}- 1 is the Doppler factor, (3 = (1 - V' 2 ) 1 / 2 , v = ZV, I'(u') is the 
specific intensity of synchrotron radiation at u', and d is the distance of the burst source. 
Here the quantities with prime are measured in the comoving frame. For simplicity we have 
ignored the relative time delay of radiation from different parts of the cone. 

For the non-relativistic blast wave, it is widely believed that the accelerated electrons 
have power-law energy distribution 71(7) oc 7~ p in the range 71 < 7 < 72, then the typical 
electron energy 71 ~ | e e^rf ^ £ /^ 2 ) where e e is the energy fraction occupied by the electrons, 
m p and m e are the mass of proton and electron. Assuming that e e is a constant, then we 
have 71 oc (3 2 . For the non-relativistic case, T ~ 1, then r = DTf3ct oc f3t/(l — /3cosO), r' = 
D(3ct oc (3t/(l — /3cosQ), the energy density u oc n/3 2 , the magnetic field strength B oc u 1 / 2 oc 
fin 1 / 2 , the peak frequency of synchrotron radiation v m = Du' m oc D^ 2 B oc / 9 5 n 1 / 2 /(l — 
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f3cosQ). Assuming that the emission spectrum I'iy') oc v'~ a , then I'iy') = V{y' )(-r-) _Q = 
I'W m ){-^)~ a oc (3 2+ba n {3+a)/2 tv~ a /{l - l3cosQ) l+a . Therefore we have the emission flux 

F(u, 9) oc u~ a (3 4+5a n i3+a)/2 t 3 g(9, (3, a) (8) 

where 

g(0, p,a)= / dtf / 3 sin9'd9'cosQ(l - pcosQ)~ (6+a) (9) 
Jo Jo 

In general, the value of g can only be calculated numerically. However here we consider the 
case 9j <C 1 and 9 <C 1, then cos© ~ cos9cos9'. In this case we can calculate the value of g 
analytically under certain conditions. After complicated calculation we find g ~ ncos9sin 2 9j, 
so finally we get 

F(u, 9) oc u- a p 4+5a n {3+a)/2 t 3 cos9sin 2 9 J (10) 



From the previous section, we can obtain the GRB afterglow light curves according to 

Q (3-3)(4+5a) + s(3+a) 

different spread velocity. For case (i), we find F(u,9) oc t 5~ s cos9sin9j, more 

clearly, F oc t ( 5 ) cos9sin 2 9j for s = 0, and F oc t~ ( r ) cos9sin 2 9j for s = 2. For case 
(ii), we find it is simple, F oc t~ ( - 5a ~ 1 ^cos9, independent of s. For case (iii), it is shown 

^ Q I 9 „ (3-s+2g)(4+5a) + 3(l-qr)(3+q) (3+2g)(5a-l) „ , , 

that F oc r + 9 5- s cos^, for s = F oc t s cos6», and for s = 2 

0,9 (l+2g)(4+5a)+2(l-g)(3+a) 

F oc t + q 3 cosfc*. It is obviously that when q — 0, it reduces to case (i), 

and when g = 1, it reduces to case (ii). 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the dynamical evolution of the non-relativistic blast wave 
in the surrounding medium (density n oc r~ s ), especially whether the sideways expansion 
will have great effect on the blast wave evolution and the GRB afterglow light curves. We 
find that whether the sideways expansion being important is strongly dependent on the 
spread velocity of the jetted material. We think that it is more reasonble for material to 
spread with sound speed, and in this case, the jet angle 9j oc lnt, so it has little effect on the 
GRB afterglow light curves, i.e. the light curves will be nearly a simple power law without a 
break, but the increase of 9j cannot be ignored, the blast wave will approach nearly spherical 
after certain time. 

In the limiting case where the spread velocity is a constant, the sideways expansion 
is very important as in the relativistic case, 9j oc t and r is nearly a constant, which is 
independent of s. In this case there should be a break in the GRB afterglow light curves. 
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However, we note that the above results is valid only when t ^> t^R, and tjvfl is usually 
about one month or even more, so we expect that the afterglow light curve should become 
steeper smoothly and gradually. As for the more general situation, we find that the light 
curve behavior should lie between the above two cases. 

Here we calculated the blast wave evolution and the jet emission features in the medium 
n oc r~ s , and especially we compare the results of s = and s = 2. It is well known that 
in the relativistic case, the light curve of s = 2 is steeper than that of s = 0, however, here 
we find that in the non-relativistic case this is not true. The slope of the light curve also 
depends on the spectral index a. For example, in case (i), we find that when a < 1.4, the 
light curve of s = 2 is steeper than that of s = 0, while if a > 1.4, the situation is opposite. 
But for case (ii), i.e. the sideways expansion is very important, we find that the light curve 
is independent of s, since in this case r ~ constant. 

Recently Frail et al. (1999) have reported on the results of 450 days observations of the 
radio afterglow of GRB 970508, and indicated that the fireball has undergone a transition to 
non-relativistic expansion at t ~ 100 days. Therefore it is very interesting and important to 
explore the dynimacal evolution of the non-relativistic jetted material in various surrounding 
medium, since we expect that over the duration of the radio emission, the emitting material 
typically become non-relativistic. 
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